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We investigate the effects of dust extinction on integrated absorption-line indices that are widely used 
to derive constraints on the ages and metallicities of composite stellar systems. Typically, absorption- 
line studies have been performed on globular clusters or elliptical galaxies, which are mostly dust-free 
systems. However, many recent studies of integrated stellar populations have focused on spiral galaxies 
which may contain significant amounts of dust. It is almost universally assumed that the effects of 
dust extinction on absorption-line measurements are entirely negligible given the narrow baseline of 
the spectral features, but no rigorous study has yet been performed to verify this conjecture. In this 
analysis, we explore the sensitivity of the standard set of Lick absorption-line indices, the higher-order 
Balmer line indices, the 4000 A break, the near-IR calcium triplet indices, and the Rose indices to 
dust absorption according to population synthesis models that incorporate a multi-component model 
for the line and continuum attenuation due to dust. The latter takes into account the finite lifetime of 
stellar birth clouds. While dust does not greatly affect the line-index measurements for single stellar 
populations, its effect can be significant for the 4000 A break or when there is a significant amount of 
current star formation. 

Subject headings: absorption-line indices, ISM: dust, extinction, galaxies: stellar content 



1. INTRODUCTION 

Stellar population studies play a fundamental role in 
our understanding of stellar evolution, initial mass func- 
tions (IMFs) associated with star cluster formation, and 
galaxy formation and evolution. With the aim of de- 
termining the star formation histories (SFHs) of stel- 
lar systems of all types (from star clusters to entire 
galaxies), many techniques have been developed, with 
varying success, to determine the luminosity-weighted 
ages and metallicities of nearby and distant stellar sys- 
tems. Important caveats hinder or limit the application 
of these techniques (e.g. Chariot, Worthey, & Bressan 
1996; MacArthur et al. 2004; Anders et al. 2004; Tantalo 
& Chiosi 2004). For example, the determination of ages 
and metallicities of stellar populations has been plagued 
by the well-known age/metallicity degeneracy (Worthey 
1994). Access to both optical and infrared imaging par- 
tially lifts this degeneracy, but broad-band colors suf- 
fer further from a degeneracy due to dust reddening 
(e.g. Bruzual, Magris, & Calvet 1988; Witt, Thronson, 
& Capuano 1992; de Jong 1996; Bell & de Jong 2000; 
MacArthur et al. 2004), which is particularly problem- 
atic in the context of gas rich, star forming systems. An- 
other method of breaking the age/metallicity degener- 
acy of stellar populations in unresolved systems involves 
measurement of surface brightness fluctuations (Worthey 
1993; Liu, Chariot & Graham 2000; Blakeslee, Vazdekis 
& Ajhar 2001), which are most sensitive to contribu- 
tions from the most luminous stars in the system - typ- 
ically evolved cool giants. Being sensitive to the sec- 
ond moment of the stellar luminosity function, surface 
brightness fluctuations provide complementary informa- 
tion to the integrated colors (the first moment of the 
stellar luminosity function). This method, however, re- 
quires the assumption of a smooth underlying light dis- 

Electronic address: |lauren@phy sics . ubc. ca| 



tribution, and thus is generally only applicable to studies 
of nearby globular clusters, ellipticals, and spiral bulges. 
This broad-band technique is also not applicable to dusty 
systems (e.g. spiral galaxies), as the resulting dumpi- 
ness will also cause fluctuations in surface brightness and 
reddening will affect the fluctuation colors. Turning to 
spectroscopy offers the possibility of overcoming the am- 
biguities of broad-band color-based analyses by studying 
the variation of individual line strengths which, if de- 
fined over a sufficiently narrow wavelength range, should 
be insensitive to the effects of dust reddening. In this 
analysis we test this conjecture by exploring the effects 
of dust absorption on commonly-used spectroscopic age 
and metallicity indicators. 

In the past three decades, tremendous progress has 
been made in the modeling of simple stellar popula- 
tions (SSPs) (e.g. Tinsley 1972; Bruzual & Chariot 1993; 
Worthey 1994; Vazdekis 1999; Fioc & Rocca-Volmerange 
1997; Maraston 1998; Bruzual & Chariot 2003). The cur- 
rent versions provide numerous observables such as high- 
resolution spectra, magnitudes, colors, mass-to-light ra- 
tios, and line-index measurements (including all those 
described above), for single bursts of star formation (SF) 
with a given IMF and metallicity at ages ranging from 
0-20 Gyr. In particular, the Bruzual & Chariot (2003) 
models include a prescription for dust attenuation (de- 
scribed in that we adopt for this analysis. 

Numerous studies of absorption-line indices in the in- 
tegrated spectra of composite systems have sought to 
disentangle and constrain the ages and metallicities of 
their stellar populations (e.g. Maraston & Thomas 2000; 
Trager et al. 2000; Schiavon et al. 2002; Caldwell, Rose, 
& Concannon 2003). The majority of these studies have 
focused on globular clusters (GCs) and elliptical galax- 
ies, where dust is not conspicuous. Recently, however, 
a number of studies have turned their attention to the 
stellar populations of spiral galaxies, where dust is an in- 
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eluctable hindrance (Fisher, Franx, & lUingworth 1996; 
Goudfrooij, Gorgas, & Jablonka 1999; Trager, Dalcan- 
ton, & Weiner 1999; Proctor & Samson 2002; Bergmann, 
J0rgensen, & HiU 2003; KaufFmann et aL 2003a; Falcon- 
Barroso et aL 2003). Reddening by dust has often been 
assumed to have a negligible effect on on line indices, but 
no rigorous study has yet been performed to verify this 
conjecture. This work presents the first such analysis. 

The outline of the paper is as follows; a description 
of the line-indices studied and their applicability as age 
and metallicity discriminators is given in 5J5] The stellar 
population synthesis models and dust prescription are 
described in SjSl In 21 we discuss the dust sensitivities of 
the Lick & higher-order Balmer line indices, the 4000 A 
break, the Ca II triplet indices, and the Rose indices, and 
the results are summarized in [jS] Finally, in Appendix IXI 
we present age and metallicity fits to the model galaxies 
in several index-index planes and discuss the errors on 
the derived physical parameters due to dust extinction. 

2. DESCRIPTION OF INDICES STUDIED 

2.1. The Lick/IDS System 

The Lick/IDS system of 21 spectral line indices was 
designed to calibrate the strength of fundamental spec- 
tral features in stars and composite systems (e.g. Gor- 
gas et al. 1993). The indices measure the strength of 
a particular spectral feature (either atomic and defined 
as an equivalent width in A, or molecular and measured 
in magnitudes) relative to a pseudo-continuum on each 
side of the feature. The most reliable indices have been 
calibrated as a function of stellar color (effective tem- 
perature), surface gravity, and metallicity (Gorgas et al. 
1993; Worthey et al. 1994) allowing for the construction 
of semi-empirical population models (e.g. Worthey 1994). 
The Lick indices are sensitive to the metallicity and age 
of stellar populations to varying degrees. When com- 
pared with population models, diagnostic plots of age 
versus metallicity sensitive indices, such as H/3 versus 
Mgb or (Fe), help break the age- metallicity degeneracy. 
However, measurements of many of the Lick indices are 
quite sensitive to spectral resolution and, thus, to the 
velocity dispersion of the system (Gonzalez 1993; Trager 
et al. 1998; Proctor & Sansom 2002), and their use re- 
quires relatively high signal-to-noise data (S/N > 50/A; 
see Cardiel et al. 1998). In addition, certain indices, H/3 
in particular, can suffer "in-filling" from nebular emis- 
sion contamination, present even in early-type galaxies 
(Gonzalez 1993; deZeeuw et al. 2002; Caldwell, Rose, & 
Concannon 2003). 

2.2. Higher Order Balmer Indices 

In order to overcome the problem of nebular emission 
fill in of the H/3 feature, Worthey & Ottaviani (1997; 
hereafter W097) introduced two pairs of indices that 
measure the higher-order Balmer lines H7 and H(5. Two 
definitions for each feature were defined, the narrower 
version (AA ~ 20 A) is denoted with the subscript "F" 
(as it encompasses all of the Balmer line absorption from 
stars of spectral type F at the Lick/IDS resolution of 8- 
10 A), e.g. H(5j7, and the wider definition (AA ~ 40 A) 
has the subscript "A" (as it includes all of the absorp- 
tion from A stars), e.g. Hja- The narrow indices are 
more age-sensitive, but require higher S/N and resolu- 



tion. While their age-sensitivity is not as strong as for 
H/3, the higher-order Balmer lines are much less affected 
by emission from ionized gases (e.g. Osterbrock 1989). 
Thus, when combined with a metallicity sensitive index, 
the W097 indices provide a more reliable age estimate 
for star-forming galaxies. Note, however, that the W097 
higher-order Balmer lines are poorly calibrated (Vazdekis 
1999), likely due to the degrading resolution at the blue 
end of the IDS data. 

2.3. The 4000 A Break 

Another widely used spectral index that obviates the 
need for high S/N and spectral resolution is the 4000 A 
break, a flux ratio that brackets the strongest discontinu- 
ity in the optical spectrum of a galaxy. The break arises 
due to the accumulation of a large number of spectral 
lines in a narrow wavelength region bluewards of 4000 
A in stellar types cooler than GO (Bruzual 1983; Gor- 
gas et al. 1999). The main contribution to the < 4000 
A opacity comes from atomic metals (Fe I, Mg I, Ca II) 
and molecular CN, which decreases for hotter and more 
metal poor stars. Thus the 4000 A break is weak for 
young and/or metal-poor stellar populations and strong 
for old, metal-rich galaxies (Kauffmann et al. 2003a). 
In its original form, (Bruzual 1983), the 4000 A break, 
denoted D(4000), was defined as the ratio of the aver- 
age fluxes per frequency unit measured over the spec- 
tral ranges 4050-4250 A and 3750-3950 A. Cardiel et al. 
(1998) demonstrated that this discontinuity can be mea- 
sured with a relative error of ^ 10% with a S/N per 
A ~ 1, thus making it better suited for lower quality 
data. However, the D(4000) does have a few drawbacks 
due to its long baseline. These were partially alleviated 
with the introduction of a narrower definition (Balogh 
et al. 1999) denoted D„(4000) and measured over the 
ranges 4000-4100 A and 3850-3950 A. The narrow defi- 
nition was designed to exploit two principal advantages: 
an improved agreement between multiple measurements 
of a given galaxy, and a weaker sensitivity to redden- 
ing effects. However, using their narrow D„(4000) index, 
Balogh et al. (1999) still had to invoke dust reddening as 
a cause for the large D„ (4000) values in a number of their 
z ^ 0.3 galaxies, i.e. D„(4000) is not impervious to dust 
effects. Nevertheless, often found in the literature is the 
statement that the D„(4000) is insensitive to dust atten- 
uation effects. For example, Kauffmann et al. (2003b) 
use the amplitude of the D„(4000) in combination with 
the strength of the US a index of W097 as diagnostics for 
the SFH of the host galaxies, from which they infer the 
dust attenuation by comparing observed to model colors, 
a method that relies heavily on the assumption of dust 
insensitivity of the D„(4000) and 116a indices. 

2.4. The Near-IR Ca II Triplet Indices 

There has been a great effort recently to extend stellar 
population studies to the near-IR region, focusing on the 
Ca II triplet as one of the most prominent features in 
the near-IR spectrum of cool stars (from spectral types 
of about F5 to M2) (Cenarro et al. 2001a,b; Vazdekis 
et al. 2003). Three "Lick-style" indices designed to mea- 
sure the strengths of the Ca II triplet lines (AA8498.02, 
8542.09, 8662.14 A), and denoted Cal, Ca2, and Ca3, 
have been defined and redefined by several authors (e.g. 
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Jones, Alloin, & Jones 1984; ArmandrofF & Zinn 1988; 
Diaz, Terlevich, & Terlevich 1989, hereafter DTT; Dehsle 
& Hardy 1992). However, measurement of a reliable con- 
tinuum for these indices is difficult due to the strong and 
crowded absorption features in the vicinity of the Ca II 
lines (largely from Fe I, Mg I, and TiO), as well as signif- 
icant blending with the hydrogen Paschen series whose 
absorption is present in stars of types G3 and hotter. 

To overcome these problems, Cenarro et al. (2001a) 
defined a new set of Ca II triplet indices specifically de- 
signed for measurement in integrated galactic spectra. 
The new definitions are categorized as "generic" indices, 
which have advantages over the classical "Lick-style" in- 
dices when looking at adjacent absorption lines in regions 
where the continuum is crowded with spectral features. 
Their CaT index includes the strenghts of all three Ca II 
lines and uses a combination of 5 continuum bandpasses 
(see Cenarro et al. 2001a for details on the measurement 
of generic indices). Additionally, they define the generic 
index PaT, which measures the strength of three of the 
H Paschen series lines that are free from Ca contamina- 
tion. Finally, the index CaT* is designed to remove the 
H Paschen line contamination from CaT making it a re- 
liable indicator of the pure Ca II triplet strength. The 
index is given by CaT* = CaT - 0.93 PaT. Vazdekis et al. 
(2003) discuss the behavior of these features for SSPs as 
predicted by recent evolutionary synthesis models. In 
the current analysis we also explore the dust sensitivity 
of the three classical Ca II triplet indices, as defined in 
DTT, and the generic indices of Cenarro et al. (2001a). 

2.5. The Rose Indices 

The optical Rose (1984, 1985) indices were developed 
in part to overcome the difficulties in identifying the con- 
tinuum in crowded spectral regions. The Rose indices 
maximize the contribution of a given feature by using 
a measurement of the ratio of its central line intensity 
to that of a close reference line, without recourse to the 
(pseudo-) continuum level. Information about equivalent 
widths is however lost in such relative measurements, 
and the Rose indices are not ideal for separating age and 
metallicity effects. Nonetheless, some of the Rose indices 
are quite sensitive to metallicity (see Vazdekis 1999) and, 
in combination with a sensitive age indicator (e.g. H/3), 
could help disentangle the two effects. The Rose indices 
also provide a unique and sensitive test for the presence 
of early-type stars (i.e. a post-starburst, and potentially 
dusty, population) in an integrated spectrum as well as 
a probe of the relative contribution of dwarf and evolved 
(red giant branch) stars to the integrated spectrum of a 
galaxy (Caldwell, Rose, & Concannon 2003). As with 
the Lick indices, it is usually assumed that the Rose in- 
dices are insensitive to dust reddening (e.g. Leonardi & 
Worthey 2000), an assumption we examine in this study. 

Finally, we also investigate dust effects on the pseudo- 
equivalent width, "Lick-style", indices of Rose (1994) and 
Jones & Worthey (1995) which have extremely narrow 
definitions (and are denoted with the subscript "HR"). 
The H7HR index has been identified as one of the most 
sensitive age distriminators (Jones & Worthey 1995), but 
these narrow features suffer a worse resolution sensitiv- 
ity than the W097 indices and thus can only be used 
with high quality data of low velocity-dispersion sys- 
tems. Vazdekis & Arimoto (1999) and Vazdekis et al. 



(2001) confronted the velocity dispersion sensitivity of 
these narrow indices by defining a set of four H7 indices 
that take resolution effects into account, allowing for re- 
liable measurements in galaxies with velocity dispersions 
up to fj ~ 300 kms^^. Note, however, that the Vazdekis 
& Arimoto (1999) indices require S/N > 200-400/ A and 
have strong error covariances due to overlapping pseudo- 
continuum and central bandpasses, making them much 
harder to measure than the W097 indices. The response 
of these indices to dust extinction is similar to the W097 
RjF and the Jones & Worthey (1995) H7HR indices, so 
the results are not shown. 

3. MODELS 

The stellar population synthesis models used for this 
analysis are those of Bruzual & Chariot (2003; here- 
after G AL AXEV) . The high resolution models using the 
Padova evolutionary tracks (Bertelli et al. 1994) and the 
Chabrier (2003) IMF were adopted. These models in- 
clude SSP spectra in the wavelength range 3200-9500 A 
with a resolution of 3 A, metallicities ranging from Z — 
0.0001-0.05 (or 0.005-2.5 times Zq), and ages ranging 
from 0-20 Gyr (in 220 unequally-spaced time steps). For 
each time step, a number of integrated quantities are 
provided such as magnitudes and colors in many differ- 
ent filter systems, line-index strengths for all definitions 
in the Lick/IDS system, W097, the two definitions of 
the 4OOOA break, and the DTT Ca II indices. The line 
strengths used here are computed directly from the high- 
resolution model spectra, i.e. the spectra have not been 
transformed to the Lick/IDS system (see e.g. W097). 
The Rose, Ca II, and "HR" spectral indices, which are 
not provided with the BC03 model outputs, were com- 
puted with a public f ortran program by A. Vazdekis^. 

The GALAXEV distribution allows for the computa- 
tion of attenuation effects due to dust according to the 
two-component model of Chariot & Fall (2000; hereafter 
CFOO). The two adjustable parameters of this model 
are fy, the total effective l^-band optical depth affect- 
ing stars younger than 10^ yr, and /i, the fraction of the 
total dust absorption contributed by diffuse interstellar 
medium (cirrus) dust. That is, 

. _ J f,/(A/5500 A)-0-^ fori<10^yr, , . 
^^~{ MTi/(A/5500 A)-0-7 for t > W yr, ^ ' 

where t is the age of any single stellar generation. This 
model has the singular feature of accounting for the fi- 
nite lifetime (~ 10'' yr) of stellar birth clouds^. The 
wavelength dependence of the effective absorption curve 
(proportional to A~"'^) was constrained to reproduce the 
observed relation between the ratio of far-infrared to ul- 
traviolet luminosities, and the ultraviolet spectral slope 
of their starburst galaxies (see also the GRASIL models 
of Silva et al. 1998). Note that t\ on the left-hand side 
of equation (QJ is a function of time and As a result, 
even if ^ = (no cirrus dust component), in the pres- 
ence of young stellar populations (t < 10^ yr), dust will 
still contribute some extinction (see Fig.^). Note that in 
the adopted dust model of CFOO, the model spectra are 



^ See http://www.iac.es/galeria/vazdckis/ 

^ Introduced to resolve the apparent conflict between the atten- 
uation of line and continuum photons in their sample of starburst 
galaxies. 
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representative of the total integrated light of a galaxy, 
i.e., neither radial nor inclination dependent information 
is available. 

Other prescriptions for dust attenuation in galaxies ex- 
ist (e.g. Witt, Thronson, & Capuano 1992; Byun, Free- 
man, & Kylafis 1994; de Jong 1996; Gordon et al. 2001), 
a few of which have been coupled with spectral syn- 
thesis and photo-ionization codes. For example, Moy 
et al. (2001) interfaced the PEGASE population syn- 
thesis models of Fioc & Rocca-Volmerange (1997) with 
the CLOUDY photo-ionization code of Ferland (2002). 
These models include a treatment for dust attenuation 
processes, but do not consider nebular emission. They 
also adopt a simple screen approximation for the dust dis- 
tribution, an unrealistic geometry for galaxies. Similarly, 
Panuzzo et al. (2003) coupled the same CLOUDY code 
with their spectrophotometric synthesis model GRASIL 
(Silva et al. 1998), to yield a complete treatment of dust 
reprocessing with more realistic geometries and account 
for the age dependence of molecular birth clouds (as in 
CFOO). While these models may be more realistic and 
provide directional information that is lacking in the 
CFOO models, the resolution of the output from these 
models is not sufficient for accurate line-index measure- 
ments. Future implementations of these models may al- 
low for a comparative analysis, but at present, the CFOO 
dust models, coupled with the BC03 SPS models, provide 
(to the best of our knowledge) the only suitable output 
for the current study. 

In a study of 705 non-Seyfert galaxies drawn from 
the Stromlo-APM redshift survey (Loveday et al. 1996), 
Chariot et al. (2002) determined parameters for the CFOO 
models that reproduce the observed integrated spectral 
properties of the nearby star-forming galaxies. These 
cover the ranges: 0.2 < Z/Z© < 4.0, 0.01 < f y < 4.0, 
and 0.2 < < 1.0, for constant and exponential SFHs 
with ages 10^ yr < i < 10^° yr. Two time-scales were 
adopted for the exponentially declining star formation 
rate: Texp = 0.1 and 6.0 Gyr. Our model realizations 
cover roughly the same parameter space. For the SSP 
models we consider values of = and 1 only since, for 
the SSP ages shown {t > 10^ yr), different values of fj, 
are equivalent to different values of fy. For the exponen- 
tial SFHs, for which there are always populations with 
t < 10^ yr present, we present models with < /i < 1 
and extend to fy = 8 (so for the fi = 0.3 case, the plots 
effectively extend to fy = 2.4). To express the effective 
dust reddening in a more familiar form, in Figure ^ we 
plot the color excesses E{B — V) and E{V — K) as a 
function of f y resulting from the dust models used in 
this analysis. The upper panels are for solar metallic- 
ity SSPs with ji = 0.0 (circles) & 1.0 (triangles) and the 
lower panels show exponential SFH with Texp = 13 Gyr 
for ^ = 0.0 (circles), 0.3 (triangles), & 0.9 (squares) at 
ages of 0.5 (dotted fines; red shades), 5 (dashed lines; 
green shades), & 13 Gyr (solid lines; blue shades). 

Before examining the response of the individual indices 
to the dust reddening, it is useful to look at the spec- 
tral energy distribution (SED) of the models. Figure |2] 
shows the model SEDs in the 3700-6700 A range for SSPs 
[left panels] and exponential SFHs with r^xp — 13 Gyr 
[right panels]. All spectra have been normalized to 5500 
A. The SSP SEDs look entirely as expected, with the 



0.5 Gyr population [top right panel] closely resembling 
an early-type (A-F) stellar spectrum with a blue contin- 
uum, strong Balmer features, but few metallic features. 
The 13 Gyr spectrum [bottom right panel] more closely 
resembles a later- type (G-K) stellar spectrum with a red- 
der continuum, weaker Balmer-line strengths, and more 
metallic absorption features. In both cases, the dust 
extincted profiles exhibit the overall low frequency red- 
dening of the SED due to the A"'''^ dependence of the 
dust model. On the other hand, the unreddened 0.5 Gyr 
SED for an exponential SFH with Texp — 13 Gyr (Fig. 12 
[top right panel, black curve]) is even bluer than its 
SSP counterpart because of the presence of young stars 
(< 0.5 Gyr) from the ongoing SF. The contrast between 
the SSP and exponential SFHs is even more striking at 
13 Gyr. The exponential SFH SED (Fig. El [bottom right 
panel, black curve]) is much bluer than its SSP counter- 
part [bottom left panel] due to the ongoing SF, but it also 
shows metallic features due to the presence of older stars 
(< 13 Gyr). In the exponential SFH case, the dust red- 
dening is more complicated due to the presence of young 
stars that still live in their birth clouds (< 10^ yr). These 
stars are more extincted than the older stars, the cirrus 
extinction component, when fi < 1.0, but being intrinsi- 
cally bright, still contribute significantly to the total flux. 
In Figure 13 [right panels] we show the case for fi — 0.3. 
Clearly, any measurement made over a long baseline (e.g. 
colors) will be affected by the dust reddening. However, 
gauging its effects on absorption-line indices from a cur- 
sory examination of the SEDs is not a straightforward 
task. In the next section we examine individual features 
and compare the dust-free indices with those computed 
with dust. 

4. RESULTS 

4.1. Dust Sensitivity of the Lick Indices 

The response of the Lick indices as a function of fy, 
/i, and age for solar metallicity {Zq — 0.02) and abun- 
dance ratio ([a/Fe] = 0) SSP models is shown in Fig- 
ure |21 In each panel, we plot Aindex versus fy, where 
Aindex is the difference between the indices measured 
with and without dust (for the same age & SFH), i.e. 
Aindex = index^Ty) — indexijv = 0). Results are 
shown for model ages of 0.5 (dotted lines; red shades), 
5 (dashed lines; green shades), & 13 Gyr (solid lines; 
blue shades). The two values for /i of 0.0 & 1.0 are de- 
noted by circles and triangles, respectively. The horizon- 
tal dotted lines in Figures |31 E & El - El denote typical 
measurement errors for the different indices; ^ 0.1 A 
for the atomic Lick indices, 0.01 mag for the molecu- 
lar Lick indices, 0.4 A for the W097 indices, 0.03 for 
D(4000) and D„(4000), and ~0.02~0.05 for the Rose in- 
dices (Gorgas et al. 1999; Jones & Worthey 1995; Kauff- 
mann et al. 2003a; Falcon-Barroso et al. 2003; Caldwell, 
Rose, & Concannon 2003; J.J. Gonzalez 2004, private 
comm.), and serve as a guide for the magnitude of the 
effect from dust required to be detectable above the noise. 
All model realizations presented here have solar metal- 
licity, Zq = 0.02. We have also considered the other 
metallicities in the GAL AXE V models. Generally, the 
dust effects are slightly more conspicuous at Z = 0.05, 
but diminish with decreasing metallicity of the stellar 
population. 

In most cases, the index measurements do not devi- 
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Fig. 1. — Color excesses E(B — V) [left panels] and E(V — K) [right panels] as a function of tv for solar metallicity SSPs with fi = 0.0 
(circles) & 1.0 (triangles) [top panels] and an exponential SFH with Texp = 13 Gyr for fi = 0.0 (circles), 0.3 (triangles), and 0.9 (squares) 
[bottom panels]. Different ages are represented by: 0.5 (dotted lines; red shades), 5 (dashed lines; green shades), & 13 Gyr (solid lines; 
blue shades). 
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ate considerably from their dust-free values in the SSP 
models, as expected. Since we consider SSPs at ages of 
0.5 Gyr and greater, by definition, models with fj, — 
have no dust extinction (thus circle symbols lie at for 
all values of fy in Fig.El). For /i = 1 models, deviations 
from the dust-free index measurements do occur. There 
are a few cases of note that show significant sensitiv- 
ity to the dust extinction, such as the molecular indices 
Mgi & Mg2, and CNi & CN2. This is not surprising 
since these indices have the longest wavelength baseline 
(AA ~400 & 200 A respectively). As evidenced by the 
molecular Mg indices, index measurements in the pres- 
ence of dust can be non-linear and age dependent, thus 
a simple "dust-correction" poses a formidable challenge. 
The effect of extinction on the Mg indices goes in the op- 
posite direction for young ages as it does for older ages 
(compare dotted and dashed/soUd Hues in Fig.|3Il. Also, 
the CN indices are more severely affected at younger ages 
whereas the Mg indices are more affected at young ages 
for small Ty S 3, but at older ages for larger fy > 4. 

Of the Lick atomic indices, those showing the great- 
est sensitivity to dust are Fe5015 & Fe4668, both 
having the longest wavelength baselines of all atomic 
indices, but their deviations from the dust-free case 
are always smaller than the typical measurement er- 
rors. In addition to the standard Lick indices, we also 
show the combination indices, (Fe)= (Fe5270-|-Fe5335)/2 
and [MgFe], first introduced by Gonzalez (1993). For 
the latter we use the slightly modified definition 
of Thomas, Maraston, & Bender (2003), [MgFe]' = 
[Mg6(0.72*Fe5270-H0.28*Fe5335)]i/2, which has the ad- 
vantage of being completely independent of the element 
abundance ratio ([a/Fe]). These combination indices 
vary very little with fy for the SSP models (Fig. O • 

The higher-order Balmer indices, shown in Figure El 
also show some sensitivity to dust at older ages in the 
SSP models, particularly the Hja index. Measurement 
errors for these indices are quite large compared to the 
other Lick indices (due to their narrow baselines), thus 
any reddening effect would likely not be detectable. 

Unlike GCs and elliptical galaxies, the assumption of 
SSP-like star formation is clearly not applicable for spi- 
ral galaxies. While the true SFHs of the latter are 
poorly known in detail, they are often approximated 
with an exponentially declining SFR, parameterized by 
a star formation timescale r^xp (e.g. Bell & de Jong 
2000; MacArthur et al. 2004). In Figure El we show 
the index responses to dust for an exponential SFH with 
Texp — 13 Gyr, which is reasonably close to a constant 
SF rate so there is a significant amount of "current" SF 
at each of the three epochs considered. 

Models with current SF yield considerably different re- 
sults than the SSPs, as seen in Figure 01 Almost all of 
the Lick indices show sensitivity to dust, some more so 
at young ages (e.g. CNi, CN2, H/3, and aU four W097 in- 
dices), others more so at older ages (e.g. Ca4227, Fe4531, 
Fe4668, Mg&, and NaD), and some show differences of 
equal magnitude, but in opposite directions at old and 
young ages (e.g. G4300, Fe4383, and the molecular Mg 
indices). The situation is undoubtedly complex and can 
be attributed to the fact that, with current SF occurring 
at all epochs, there is a significant and non-linear amount 
of extinction of the hottest, youngest stars which con- 



tribute considerably to the total optical flux. As such, 
dust effects are often stronger for models with smaller fj, 
and these effects saturate at optical depths above about 
f y — 3 (at which the birth clouds become completely 
optically thick). 

The situation for mixed stellar populations with signif- 
icant amounts of current SF is clearly complicated and 
not easily prescribed. However, for many of the indices 
(e.g. Ca4227, Ca4455, Fe5406, Fe5709, Fe5782, NaD, 
TiOi, and Ti02), the magnitude of the dust effects is 
never close to that of current typical measurement er- 
rors, and thus would not be detected. Unfortunately, 
however, the most often used indices (e.g. II/3, Mg6, the 
composite indices (Fe) and [MgFe]', and the W097 in- 
dices) are also the most seriously affected. 

It is useful to take a closer look at the cause for the 
erratic behavior in some of the indices in response to the 
dust extinction. In Figure we focus on the II/3 index. 
Plotted in each panel are solar metallicity SEDs for fy 
— (black), 3 (gray), & 6 (pink) at ages of 0.5 [top 
panels], & 13 Gyr [bottom panels]. The index red and 
blue pseudo-continuum band limits are marked by the 
green dashed vertical lines and the central band by the 
solid vertical lines. The dotted lines mark the pseudo- 
continuum level for each spectrum. The spectra are 
all normalized to their 4861 A flux (H/? line center). Fig- 
ure plots SSPs [right panels] with /i = 1.0 exponential 
SFHs [left panels] with /z = 0.3. There is very little ef- 
fect in the SSP index measurement, the change in slope 
of the SED is matched by the change in slope of the 
pseudo-continuum. However, for the exponential SFH at 
0.5 Gyr (Fig. [top right panel]), the continuum is in- 
creased relative to the depth of the absorption-line in the 
dust extincted models. This can be understood by rec- 
ognizing that the young OB stars that suffer significant 
extinction have lower H/3 values. Thus, by hiding these 
low H/3 index stars, the index is effectively increased. In 
fact, this is the case for all of the Balmer lines, hence all 
of the Balmer-line indices show similar behavior. This 
also highlights a major limitation when using Balmer- 
line indices in determining young ages (< 0.5 Gyr). The 
models become degenerate in age here because there are 
two possible ages that can be inferred at these Balmer 
index values, i.e. the Balmer indices are double- valued in 
the region of young stellar population ages. This point 
and its resulting limitations are further emphasized in 
Appendix 1X1 

For the molecular indices, which have much broader 
baselines and crowding of spectral lines in their passband 
regions, the situation can be more complicated. For in- 
stance, in Figure El we focus on the Mgi index. The 
color codes, line-types, and SFHs are the same as in Fig- 
ure |5l In this case, given the much broader baseline, the 
slope of the pseudo-continuum in the extincted models 
changes quite dramatically. In the 0.5 Gyr SSP (Fig. El 
[top left panel]) and 13 Gyr exponential SFH (Fig. El 
[bottom right panel]) cases it changes from being a nega- 
tive slope in the dust-free model (black) to positive in the 
most extincted curve (pink). These pseudo-continua pass 
through a crowded set of absorption lines in the central 
passband, and a simple inspection of these figures does 
not reveal which model would yield larger index measure- 
ments. It is clear from these figures that a simple dust 
correction prescription for all of the indices is cannot be 
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Fig. 4. — Same as Figure l3l but for for an exponential SFH with Texp = 13 Gyr and values of are represented as 0.0 (circles), 0.3 
(triangles), and 0.9 (squares). 
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Fig. 5. — Comparison of the H/3 index for solar metallicity SSP models with fi = 1.0 [left panels], and exponential SFH Texp = 13 Gyr 
models with fi = 0.3. [left panels]. All panels show models with ry = (black), 3 (gray), & 6 (pink) at ages of 0.5 [top panels] & 13 Gyr 
[bottom panels]. The index red and blue pseudo-continuum band limits are marked by the green dashed vertical lines and the central band 
hy the solid vertical lines. The dotted lines mark the pseudo-continuum level for each spectrum. The spectra are all normalized to their 
flux at 4861 A (the center of the H/3 line). 
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Fig. 6. — Same as Figure O but for Mgi. The spectra are all normalized to their flux at 5107 A (roughly the center of the index central 
passband) . 
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achieved. 

4.2. Dust sensitivity of the 4OOOA Break 

Figure 13 shows the response of D(4000) [top panels] 
and D„(4GGG) [bottom panels] to attenuation from dust 
for SSPs [left panels] and exponential SFHs with r^xp = 
13 Gyr [right panels]. The narrower index is certainly less 
affected by dust (note the different y-axis scales), but it 
is by no means impervious to dust reddening effects, and 
can reach deviations with magnitudes much larger than 
the typical error for this index. Dust effects are smaller 
at younger ages, but they can be significant at all ages 
for realistic values of fy. A detailed discussion of the 
influence of these dust effects on the determination of 
ages and metallicities is deferred to appendix IXI 

Gorgas et al. (1999) suggest that D(4G0G) could be 
corrected for the effects of internal reddening according 
to, 

D(4000)_,,, = D(4000)„,_, * 10^o.0988i^(B-y)^ 

(2) 

derived using the mean extinction curve from Savage & 
Mathis (1979). However, this assumes that the dust is 
distributed as a screen, which is likely unrealistic for 
galaxies (e.g. Witt, Thronson, & Capuano 1992), and 
that the value of the color excess E{B — V) is known. 
Even knowing E{B — V), as we do for these models, the 
above formula does not reproduce the dust-free D(4000) 
value for the GFGG dust model. The numerical constant 
in equation ^ could be adjusted to obtain a good fit, 
but this number is a function of age, depends on the dust 
properties, and requires a priori knowledge of E(B — V). 
Hence, we do not explore this putative correction further, 
but note that, as we concluded for the Lick indices, there 
is likely no such simple dust correction for the 4000 A 
break of mixed stellar populations and complex dust ge- 
ometries. 

4.3. Dust Sensitivity of the Ca 11 Triplet Indices 

The response of the Ca II triplet indices of DTT and 
Cenarro et al. (2GGla) as a function of fy, /U, and age for 
solar metallicity SSP models are shown in Figure |H1 and 
for an exponential SFH with r^xp = 13 Gyr in Figure El 
The three classical Lick-style Ca indices of DTT (left) 
show some sensitivity to dust extinction for the SSPs, 
Ca3 showing the largest sensitivity, but the deviations 
remain below the typical measurement errors thus would 
not be detected. The generic CaT indices (right) are vir- 
tually unaffected by dust extinction in the SSP models 
at any age. As usual, the situation is not as straightfor- 
ward for the exponential SFH (Fig. El . Here, the only 
index that is not significantly affected is the Gal index of 
DTT. All of the generic indices of Cenarro et al. (2GGla) 
are significantly altered due to dust extinction, particu- 
larly at young ages. As in the case for the Lick indices in 
the exponential SFH models f ii4.ll Fig. 0, dust effects 
are often stronger for models with smaller /.t and the ef- 
fects saturate at optical depths above about fy = 4. 
In general, though, for models that do not have signif- 
icant amounts of current SF, the dust reddening effects 
in the Ca II triplet indices are still small compared to 
measurement errors, and should not cause problems in 
their interpretation. 



4.4. Dust Sensitivity of the Rose Indices 

The response of the Rose (1984, 1994) and Jones 
& Worthey (1995) indices are shown in Figure ^] for 
SSPs and in Figure ^2 for an exponential SFH with 
Texp — 13 Gyr. Symbols and line types are as in Fig- 
ures 13 & 0] For SSP models there are a few indices that 
are considerably modified by dust reddening, namely 
the H(5/Fe IA4045, H5/Fe IA4G63, Sr II/Fe IA4G45, 
H7/Gband, and Ca II indices, particularly at older ages. 
All other indices, including the pseudo-equivalent width 
indices, are essentially unaffected by the dust reddening 
in the SSPs at any age. For the exponential SFH, shown 
in Figure [TTl the response of the indices gets quite com- 
plicated. The offset from the dust-free case can be either 
positive (typically for the older models) or negative (typ- 
ically for youngest models). Again, for the fi = G.O mod- 
els and for some of the ^ — 0.3 models, the index offsets 
saturate for optical depths above f\/ = 3. Ultimately 
though, even with current SF, for the majority of the 
Rose indices, the index offsets due to dust reddening are 
modest. On the other hand, with a significant amount of 
current SF, the pseudo-equivalent width indices can be 
quite affected, more so at young ages for H7HR, but more 
so for old ages for Ca Ihr and Fe Ihr- The "HR" indices 
are virtually unaffected for SSP models at all ages. 

5. DISCUSSION 

As an illustration of the potential dangers of dust ex- 
tinction when using absorption-line indices to determine 
ages and metallicities of star forming galaxies, in Fig- 
ure El we present model tracks for the Lick H/3 versus 
(Fe) diagnostic plot for the dust-free exponential SFH 
with Texp = 13 Gyr for Z = G.0004, G.02, and 0.G5 (gray 
curves). The circles are the index values for the same 
SFH, but with dust extinctions ranging from fy = G 
(largest point size) to f y = 8 (smallest point size). There 
is virtually no error for old and metal-poor galaxies, but 
significant problems can arise at higher metallicities. For 
example, the location in the H/3 versus (Fe) plane for a so- 
lar metallicity stellar population at any age with f y — 1 
lies closer to the Z = G.05 model curve. The age er- 
rors are not significant in this case, but as extinction 
increases, the data points lie significantly outside the re- 
gion covered by the model grids. 

However, since in reality we don't know the true SFH 
of a galaxy, what is often computed are "SSP-ages and 
metallicities" , i.e., the values that would be measured if a 
given index-index combination is fit to SSP model grids. 
In order to gauge the magnitude of the errors on SSP 
ages and metallicities due to dust extinction, we fit the 
model galaxies to SSP grids of different index-index com- 
binations. The error on the derived parameters is taken 
as the difference between the fit for the dusty model and 
the corresponding dust-free model (i.e. AA = Age(fy) 
— Age(f y = G)). In some cases we encounter difficulties 
with the models with current SF, as these models can 
lie off the SSP model grids, even in the dust-free mod- 
els. Additionally, due to the double-valued nature of the 
Balmer lines at young ages (discussed in M4.1(l , there is an 
additional degeneracy for points that lie in the youngest 
region of the SSP model grids. This was the case for 
all of the exponential SFH models at G.5 Gyr, for which 
a large fraction of the stars are very young. Details of 
the fits and tables providing the resulting age and metal- 
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Fig. 7. — D(4000) [top panels] and Dn{4000) [bottom panels] differences, Aindex = indexirv) — index{Ty = 0), as a function of ty 
for SSPs [left panels] and for an exponential SFH with r^xp = 13 Gyr [right panels]. Different ages are represented by: 0.5 (dotted lines; 
red shades), 5 (dashed lines; green shades), &: 13 Gyr (solid lines; blue shades). The black horizontal dotted lines represent the typical 
measurement error. 
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Fig. 8. — Calcium triplet index differences, Aindex = 
indexirv) — index{Tv = 0), as a function of tv for solar metallic- 
ity SSPs. Different ages are represented by: 0.5 (dotted lines; red 
shades), 5 (dashed lines; green shades), & 13 Gyr (solid lines; blue 
shades). Different values of fi are represented as 0.0 (circles), and 
1.0 (triangles). The black horizontal dotted lines represent typical 
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Fig. 9. — Same as Figure ISl but for for an exponential SFH 
with Texp = 13 Gyr. Different values of fi are represented by 0.0 
(circles), 0.3 (triangles), and 0.9 (squares). 
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licity errors for a selected number of index-index diag- 
nostic plots are provided in Appendix ^ The primary 
result from this analysis is that, in general, when dust 
extinction affects the SSP age and metallicity determi- 
nations, the errors on the physical parameters are of the 
same order as the la confidence intervals from typical 
measurement errors, and thus would not likely be de- 
tected above the noise. The most notable exception is 
any index-index plane using the D„(4000), which suffers 
significant errors even in the SSP models. Index-index 
grids using the longest baseline indices can also suffer 
significant errors on the measured physical parameters 
due to the effects of dust. 

Figure [T51 illustrates the potential pitfalls of discount- 
ing dust extinction in an analysis that uses index di- 
agnostics to determine SFHs. There, we plot H5a ver- 
sus Dji(4000) for the dust- free SSP (gray curves) and 
exponential SFH (purple curves) models for solar and 
Z = 0.05 metallicity and ages up to 20 Gyr. The tri- 
angles represent the /i = 0.3 model indices for fy = 
(largest point size and darkest shade) to f\/ = 8 (smallest 
point size, lightest shade). The locus of galaxies in the 
H^A - D„(4000) plane has been used to determine their 
SFHs (e.g. Kauffmann et al. 2003a). Clearly, if a galaxy 
contains a significant amount of dust, its location in this 



plane will be best fit by a SFH that is different from its 
true SFH. As shown, the exponential SFHs move closer 
to the SSP region when dust is added. 

Fortunately, in the Rose dwarf/giant diagnostic plot 
(Sr II A4077/Fe IA4063 versus H5/A4063), the dust ef- 
fects run along the lines of constant surface gravity. 
Thus, this diagnostic plot is still useful for determening 
the relative amounts of young versus evolved stars, even 
when dust is present. 

It must be emphasized that our analysis has not ac- 
counted for emission-line contamination. Dust effects 
aside, this could potentially be the most worrisome pol- 
lutant for all age determinations involving Balmer lines, 
particularly for the youngest stellar populations. While 
considering the higher-order Balmer lines (H7 and H(5) 
greatly reduces the problem, it does not go away alto- 
gether. 

To summarize, we have considered dust effects on 
absorption-line indices using the Bruzual & Chariot 
(2003) population synthesis models incorporating the 
multi-component model of CFOO for the line and con- 
tinuum attenuation due to dust. For quiescent stellar 
populations (e.g. spheroids and globular clusters), dust 
extinction effects are small for most indices. A notable 
exception is the 4000 A break, whose sensitivity to dust 
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13 Gyr. Different values of fj, are represented by 0.0 (circles), 0.3 



can translate into significant errors in the age determi- 
nation of the stellar population. For models with cur- 
rent SF, many of the indices are significantly modified 
due to dust reddening effects, and their behavior de- 
pends on age, dust distribution, and the effective opti- 
cal depth. Unfortunately, no simple dust-correction can 
be prescribed, but future spectroscopic studies of stellar 
populations in dusty environments (e.g. late-type spirals) 
ought to consider possible effects due to dust attenuation 
in their measurements and its resulting effect on physical 
interpretations. Fortunately, there are particular indices 
that are negligibly affected by dust extinction (compared 
to the current level of measurement precision) at different 
regions in the physical parameter space, and the safest 



current approach consists of mapping a broad range of 
indices with the shortest wavelength baselines. 
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curves), Z = 0.0004 (triangles), 0.02 (circles), and 0.05 (squares). 
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10.0, and 20.0 Gyr. The colored points (red, green, and blue shades 
for 0.5, 5, and 13 Gyr respectively) are the model indices measured 
for the same exponential SFH models, but with dust extinctions 
from TV = (largest point size) to ry = 8 (smallest point size). 
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APPENDIX 

AGE AND METALLICITY FITTING 

In order to translate absorption-line indices to a physical age and metallicity scale, the index measurements must 
be compared with stellar population synthesis models. In their most basic form, commonly referred to as simple 
stellar populations (SSPs), these models provide evolutionary information for a coeval population of stars born with 
a given composition and initial mass function (IMF). Several such SSP models have been produced by a number of 
independent groups and are in a constant state of flux as improvements to many of the input parameters (e.g. stellar 
libraries, model atmospheres, convection, mass loss, mixing) come to light. There are discrepancies among the different 
models that, depending on the application, may result in significantly different interpretations of the observed stellar 
population signatures. We are not interested here in a detailed comparison of the different SPS models, but rather 
seek a quantitative guide for the magnitude of the errors due to dust effects on the physical parameters derived from 
absorption-line indices. To this end, we use the 2003 implementation of the Bruzual & Chariot (2003) SPS models 
(GALAEX) for the age and metallicity determinations, the same as are used to compute the model "galaxies" (see 

®- . . . 

While most galaxies have likely undergone SFHs that are much more complex than the simple SSP described above, 

the common practice for line-index measurements is to compute "SSP-ages and metallicities" , that is, compare the 
galaxy indices to the SSP grids, irrespective of the true SFH of the galaxy (which is unknown). We adopt this technique 
here and thus refer to SSP parameters even for the models with exponential SFHs. 

Ages and metallicities are determined by fitting the model galaxy index measurements to the SSP model grids using 
a maximum-likelihood approach. We retain the same grid in time as is provided in the model SSPs (220 unequally 
spaced time steps from t = 0-20 Gyr) and interpolate (linearly) between the model metallicities on a fine grid of 
120 unequally spaced steps in Z. In order to compare all indices on a similar scale, the atomic indices computed as 
equivalent widths in A, Iew, are converted to a magnitude, Imag, as: 

I^aa = -2.5 logio (l - ^) (Al) 

where AA is the width in A of the index feature bandpass. The corresponding error in magnitude is: 

2.5 1ogio(e) 

^ AA lO-O-^^—a ^ ' 

The two definitions of the 4000 A break, which are dimensionless flux ratios, are converted to a magnitude scale as: 

D4000„,a<, = -2.5 logio (mOOOraUo) (A3) 
with the corresponding error in magnitude as: 

,CT[D4000™t^] 
D4000™t»o 

We compute an age and metallicity by minimizing the following figure of merit: 

i—1 



a[D4000_,] = 2.5 logio(e) Trn" ■ (A4) 



where N is the number of indices (only 2 are used here, but this method could be generalized to fit multiple indices 
simultaneously), Oi is the "observed" value (in magnitudes) of index i and SOi is its error in magnitudes, and M{A, Z\ 
is the SSP model value of index i for a given age and metallicity combination. 

Errors for the individual age and metallicity measurements are estimated using a Monte Carlo approach. One 
thousand realizations of the model fits are computed using errors drawn from a normal distribution of the observational 
errors (taken here as the "typical" observation errors for each index quoted in ti4.1|) . The errors for the measured ages 
and metallicities are taken as half the interval containing 68% of t he 1000 Monte Carlo realizations (i.e. the Icr 
confidence interval). Two examples of the fits are given in Figure lA 1 41 which shows the logarithmic difference between 
the physical parameters derived from each simulated point r elativ e to the best- fit value for the D„(4000) - Fe4668 
plane (see also the corresponding index-index plot in Figure IAl5l [right panel]). The left panel presents the fit for 
the dust-free (/i = 0.0) SSP case which indeed finds the correct model age (A = 5 Gyr; Alog]^Q[Age (Gyr)] = 0) and 
metallicity (Z = 0.02; A[log]^Q(Z/Z0)] ~ 0), with corresponding Icr confidence intervals (represented as error bars) of 
~ 0.2 dex for both parameters. The right panel presents the fit for the dusty case (/i = 1.0) at 0.5 Gyr (but note the 
different y-axis scales). The f y = fit finds the correct physical parameters, but as the dust opacity increases, the 
best fit wanders around in age (from A — 0.5 Gyr at f y = to A = 10.0 Gyr at Ty — 6) and metallicity (from Z = 
0.02 at fy = to Z = 0.003 at fy ~ 6). Since we do not consider extrapolations, there are "edges" in these plots 
that represent the limits of the model grids. As a result, the confidence limits can be underestimated for points fit 
near the model grid limits, which also generally correspond to a region of greater degeneracy in age and metallicity. 
As a result, any points that lie near the model limits should be treated with caution. 
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Given that we are comparing model galaxies with exponential SFHs against SSP grids, it is inevitable that some of 
the parameter space will fall outside the model grids and thus cannot be fit reliably. In fact, even for some of the SSP 
models, the dust extinction can be severe enough to push the locus of the model point in the index-index diagram 
off the SSP model grids. Both of these situations are observed in the D„(4000) versus Fe4668 index-index diagram in 
Figure fAl5l [right panel]. 

Tables ^ & 121 present the results from our age and metallicity fitting for a selected number of index-index combi- 
nations: H/3 versus Mg6, Mg2, and [MgFe]' (top) and H/3, RSa, and D„(4000) versus Fe4668 (bottom) for the SSP 
and exponential SFH model galaxies respectively (tables for other index combinations are available from the author 
upon request). For each index-pair, we list the model age, MA, the effective F-band optical depth, Ty, whether the 
fit is reliable in the "fit" column (see below), the best-fit age with its la error in brackets, the difference between the 
measured age in the extincted models and the dust-free age, AA = Age(fy) — Age(fy = 0), and the same for the 
metallicities, which are given as log2o(^/-^0) ^^'^ = ^^Sioi^ / ^Q)i'^v) — logio(^/Z0)(fy = 0). Reliable fits are 
indicated with a check mark, ■/ . For unreliable fits the check mark is replaced with a code indicating the reason for 
which they could not be fit. The codes are: "OG" if the point lies off the model grids, "DB" if it lies in the young 
double-valued Balmer line region, and "DM" if it lies in a degenerate metallicity region (only occurred for the young 
exponential SFH model in the D„(4000) - Fe4668 plane). Note that, in the absence of a reliable fit for the dust-free 
case, even if reliable fits exist for the corresponding dusty models, the A's are not well defined (only occurred for the 
0.5 Gyr exponential SFH model in Tabled in the D„(4000) - Fe466 8 plan e). Examples of the messy double- valued 
Balmer line region can be seen in the H/? versus Mg6 plot in Figure fA 151 [left panel] where the ages plotted extend 
down to 0.002 Gyr, as weU as in Figures [HI fc IIH in fgl 

It is obvious from Tables (and see Fig. lA14ll that there is significant uncertainty in the fits given our assumed 

"typical" measurement errors, regardless of the dust content. For example, the 0.5 Gyr SSP model la confidence 
limits (Table nj range from '^0.1-0.4 Gyr in age and ~0.2-0.3 dex in metallicity, depending on the index-index grid 
used. Note that the la confidence limits are not only dependent on the particular index-index plane used, but they 
also depend on the location of the point within the grid (in other words, to what degree the age and metallicity 
are separated at the given locus). In general, the degeneracies become more severe at the extreme metallicities (i.e. 
Z < 0.001 and Z > 0.02), and at young ages, where both the narrower region covered by the metallicity span and the 
double-valued Balmer lines contribute to the problem. 

As expected from 21 the fits for the 0.5 Gyr SSP models (Table^l are virtually unaffected by dust (for the index- 
index planes shown), with the exception of the D„(4000) - Fe4668 plane, for which the age and metallicity errors are 
large (up to 9.5 Gyr and 0.8 dex respectively) and significant compared to the la confidence limits, even for small 
amounts of dust. Significant errors occur at the highest dust opacities for the 5 Gyr models (in all but the H/3 - Mgfe 
and [MgFe]' planes), but in this case they are of the same order as the la confidence limits. The errors due to dust on 
age and metallicity are most severe in the 13 Gyr models, with some points ending up off the model grids at the highest 
fv's (e.g. in the H/3 and USa versus Fe4668 grids), but these are again matched by the large la confidence limits for 
fits in these regions of the grid. In summary, for the SSP model galaxies, if the model indices lie in a region of the 
index-index diagnostic plot where the age and metallicity are well separated (close to orthogonal) , the fits are generally 
independant of dust extinction for most index combinations. Where the fits are affected by the dust (generally at 
older ages and for indices with broader baselines), the errors in the derived ages and metallicities are generally of the 
same order as the la confidence limits and would thus be difficult to detect. The most notable exceptions here are the 
H/3-Mg2 plane, where the dust pushes the points for the 0.5 Gyr models into the degenerate double-valued Balmer 
line region, and the D„(4000) - Fe4668 plane, where the errors due to dust become large at young ages, and at old 
ages the points extend off the edges of the model grids. 

On the other hand, for the exponential SFH models, again as expected from f^l the situation is more complicated. 
Examination of Table |21 reveals that none of the 0.5 Gyr models can be fit reliably in the Balmcr-index versus 
metallicity-index grids as they all fall in the double- valued Balmer line region of the plots. Reliable fits for the dustier 
0.5 Gyr models in the D„(4000) - Fe4668 plane arc obtained, but here the A's are undefined due to the poor fit in the 
dust-free model. Reliable fits are found for most of the 5 Gyr exponential SFH models with "SSP ages" in the range 
~0. 9-1.0 Gyr (i.e. this is the luminosity-weighted SSP age of a 5 Gyr old stellar population that has been forming 
stars at a constant rate, roughly, throughout its lifetime). The ages, when reliably fit, generally agree between the 
different index-index plane fits. When dust has an effect on the ages, it tends to make them slightly youger, but the 
AA's are always of the same order as the la confidence limits for the fits. The metallicity fits for these models, which 
are all inherently solar metallicity, are in the range logiQ{Z/ZQ) ^ —0.7 to -1-0.2 (or Z = 0.004 - 0.03). They are 
generally lower(higher) than solar in the 5 (13) Gyr models and increase with dust extinction. Unlike the age fits, the 
metallicity fits do not agree very well between the different index-index planes. For example, the fy = fit in 13 Gyr 
models gives \ogiQ{Z/ZQ) = -1-0.32 and —0.12 dex in the H/3 - Mgb and H/3 - Fe4668 planes respectively. The ages for 
the 13 Gyr exponential SFH model fits are also in the '^0.9-1.0 Gyr range. While these do have weaker Balmer lines 
than the 5 Gyr models, their metallicity index values are slightly stronger. This has the overall effect of moving them 
along a line of roughly constant age. An example of this is shown by the 5 Gyr (green) and 13 Gyr (blue) triangles in 
the H/3 versus Mg6 diagnostic plot in Figure 1^151 [left panel]. These results highlight once again the greater challenge 
of fitting reliable ages and metallicities for "young" populations. 
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Fig. A14. — Examples of the Monte Carlo fits for D„(4000) versus Fe4668 diagnostic plots, represented as the logarithmic difference 
between the p hysic al parameter derived from each simulated point relative to the best-fit value (also see the corresponding index-index 
plot in Figure IA15I right panel). The error bars represent the la confidence limits. Left: SSP model at 5 Gyr with fi = 0.0. Right: SSP 
model at 0.5 Gyr with fi = 1.0 (note the different y-axis scales). 




Fig. A15. — H/3 versus Mg6 [left] and D„(4000) versus Fe4668 [right] diagnostic plots. The grids are dust-free BC03 SSP models. Lines 
of constant age (dashed) are shown for ages 0.2, 0.5, 1.0, 2.0, 5.0, 10.0, and 20.0 Gyr in both panels, but in the H/3 versus Mg6 plot (left), 
the age extends to 0.002 Gyr to illustrate the messy young age, double-valued Balmer line, region (see text). Lines of constant metallicity 
(solid) are shown for Z = 0.0001 (purple), 0.0004 (dark blue), 0.004 (light blue), 0.008 (green), 0.02 (yellow), and 0.05 (red). The circles 
are the indices measured from the solar metallicity SSP models with /i = 1.0 for Ty = (largest point size) to fy = 6 (smallest point 
size), and the triangles are the exponential SFH model indices with Texp = 13 Gyr and fi = 0.3 for fy = (largest point size) to fy = 8 
(smallest point size). The red, green, and blue shades are for model ages of 0.5, 5, and 13 Gyr respectively. 
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1 in brackets next to the age and metallicity measurements are the Icr confidence limits 

distribution of the nicasurcnicnl errors. 



Note. — All ages are given in Gyr and metallicities as \o^-^q{Z / Zq). The numbers i 

on the derived values based on 1000 Monte Carlo realizations drawn from a Gaussian 

^A check mark, •/ , indicates a reliable fit. Unreliable fits are coded as: "OG" if the point lies off the model grids, and "db" if it lies in the young double- valued Balmer 
line region. 
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TABLE 2 

Age iz Metallicity Errors: exponential SFH with rexp = 13 Gyr and // = 0.3. 
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0.60) 0.58 -0.60 

0.27) 0.02 0.02 

0.54) -0.03 0.40 

0.26) 0.03 -0.05 



0.52) 0.00 

0.53) -0.51 

0.53) 0.12 

1.08) 0.49 

0.51) 0.05 



0.06 
0.64 
0.57 
0.57 
0.26 



0.,33) 0.00 

0.54) 0.58 

0.19) 0.51 

0.19) 0.51 

0.39) 0.20 



-0.01 
-0.82 
-0.76 
-0.76 
-0.70 



0.43) 0.00 

0.19) -0.82 

0.10) -0.75 

0.08) -0.75 

0.09) -0.69 



1.02 
1.02 
0.81 
0.81 
0.81 



0.11) 
0.21) 
0.21) 
0.21) 
0.21) 



0.00 -0.18 

0.00 -0.68 

-0.21 0.20 

-0.21 0.20 

-0.21 0.20 



0.45) 0.00 

0.90) -0.50 

0.41) 0.38 

0.36) 0.38 

0.34) 0.38 



0.90 
0.81 
0.81 
0.81 
0.81 



0.10) 0.00 0.10 

0.10) -0.10 0.22 

0.10) -0.10 0.20 

0.09) -0.10 0.20 

0.10) -0.10 0.20 



0.19) 0.00 

0.12) 0.12 

0.12) 0.11 

0.12) 0.11 

0.12) 0.11 



1.02 
0.90 
0.90 

0.90 
0.90 



0.00) 0.00 -0.22 

0.11) -0.11 0.04 

0.10) -0.11 0.04 

0.10) -0.11 0.04 

0.10) -0.11 0.04 



0.19) 
0.22) 
0.22) 
0.21) 
0.22) 



0.90 
0.90 
0.90 
0.90 
0.90 



0.11) 0.00 0.32 

0.10) 0.00 0.33 

0.10) 0.00 0.33 

0.10) 0.00 0.33 

0.10) 0.00 0.33 



0.07) 0.00 

0.04) 0.01 

0.04) 0.01 

0.04) 0.01 

0.04) 0.01 



0.90 
0.90 
0.90 
0.90 
0.90 



0.11) 0.00 0.33 

0.10) 0.00 0.33 

0.10) 0.00 0.33 

0.10) 0.00 0.33 

0.10) 0.00 0.33 



0.08) 0.00 

0.04) 0.00 

0.04) 0.00 

0.04) 0.00 

0.04) 0.00 



1.02 
0.90 
0.90 
0.90 
0.90 



0.12) 0.00 0.20 

0.11) -0.11 0.32 

0.10) -0.11 0.32 

0.10) -0.11 0.32 

0.10) -0.11 0.32 



0.18) 0.00 

0.07) 0.12 

0.07) 0.12 

0.06) 0.12 

0.06) 0.12 



vs. Fe4668 



H5a vs. Fe4668 



D„(4000) vs. Fe4668 



DB 
DB 
DB 
DB 
DB 



0.04 

0.64 
0.64 
0.64 
0.64 



0.34) 0.00 0.22 

0.58) 0.60 -0.60 

0.57) 0.60 -0.54 

0.57) 0.60 -0.54 

0.57) 0.60 -0.54 



0.44) 0.00 

0.72) -0.82 

0.33) -0.76 

0.32) -0.76 

0.32) -0.76 



DB 
DB 



0.04 
0.81 
0.81 
0.81 
0.81 



0.98) 0.00 0.28 

0.42) 0.77 -0.72 

0.72) 0.77 -0.72 

0.72) 0.77 -0.72 

0.72) 0.77 -0.72 



0.99) 0.00 

0.76) -1.00 

0.74) -1.00 

0.74) -1.00 

0.72) -1.00 



DM 
OG 
■/ 



0.07 
0.18 
0.20 
0.26 
0.32 



0.05) 0.00 -0.01 

0.11) 0.11 0.40 

0.04) 0.13 0.37 

0.05) 0.18 0.10 

0.07) 0.25 -0.22 



0.32) 0.00 

0.80) 0.41 

0.32) 0.38 

0.39) 0.11 

0.34) -0.21 



1.02 
1.02 
1.02 
1.02 
1.02 



0.00) 0.00 -0.24 

0.11) 0.00 -0.24 

0.11) 0.00 -0.24 

0.11) 0.00 -0.24 

0.11) 0.00 -0.24 



0.07) 0.00 

0.12) 0.00 

0.13) 0.00 

0.14) 0.00 

0.12) 0.00 



1.02 
0.90 
0.90 
0.90 
0.90 



0.12) 0.00 -0.25 

0.10) -0.11 -0.09 

0.15) -0.11 -0.09 

0.15) -0.11 -0.09 

0.15) -0.11 -0.09 



0.14) 0.00 

0.16) 0.17 

0.24) 0.17 

0.23) 0.17 

0.24) 0.17 



OG 
OG 
OG 
OG 
OG 



0.32 
0.45 
0.45 
0.51 
0.51 



0.03) 0.00 0.40 

0.05) 0.13 0.39 

0.05) 0.13 0.40 

0.06) 0.19 0.38 

0.06) 0.19 0.40 



0.00) 0.00 

0.01) -0.01 

0.00) 0.00 

0.02) -0.02 

0.05) 0.00 



1.14 
1.02 
1.02 
1.02 
1.02 



0.12) 0.00 -0.12 

0.12) -0.12 0.18 

0.12) -0.12 0.18 

0.12) -0.12 0.18 

0.12) -0.12 0.18 



0.20) 0.00 

0.20) 0.30 

0.19) 0.30 

0.19) 0.30 

0.19) 0.30 



OG 
OG 
OG 
OG 



1.02 

0.81 
0.81 
0.81 
0.81 



0.17) 0.00 0.06 

0.33) -0.21 0.37 

0.33) -0.21 0.37 

0.33) -0.21 0.37 

0.33) -0.21 0.37 



0.26) 0.00 

0.22) 0.31 

0.21) 0.31 

0.21) 0.31 

0.21) 0.31 



OG 
OG 
OG 
OG 
OG 



0.57 

1.02 
1.02 
1.02 
1.28 



0.07) 0.00 0.40 

0.37) 0.44 0.00 

0.28) 0.44 0.00 

0.26) 0.44 0.00 

0.26) 0.71 0.00 



0.00) 0.00 

0.40) -0.40 

0.34) -0.40 

0.00) -0.40 

0.02) -0.40 



Note. — All ages are given in Gyr and mctallicitics as log^Q(.Z/ZQ). The numbers in brackets next to the age and metallicity measurements are the \a confidence 
limits on the derived values based on 1000 Monte Carlo realizations drawn from a Gaussian distribution of the measurement errors. Note that the A's are not well 
defined for the 0.5 Gyr models in the Dn(4000) vs. Fe4668 that are marked as reliable fits, since the dust-free case cannot be reliably fit. 

^ A check mark, / . indicates a reliable fit. Unreliable fits arc coded as: "OG" if the point lies off' the model grids, ''db" if it lies in the young double- valued Balmer 
line region, and "dm" if it lies in a degenerate metallicity region (only occurred for the young exponential SFH model in the D„(4000) — Fe4668 plane). 



